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We develop a mean-field theory for random quantum spin systems using the spin coherent state 
path integral representation. After the model is reduced to the mean field one-body Hamiltonian, 
the integral is analyzed with the aid of several methods such as the semiclassical method and 
the gauge transformation. As an application we consider the Sherrington-Kirkpatrick model in a 
transverse field. Using the Landau expansion and its improved versions, we give a detailed analysis 
of the imaginary-time dependence of the order parameters. Integrating out the quantum part of 
the order parameters, we obtain the effective renormalized free energy written in terms of the 
^ ' classically defined order parameters. Our method allows us to obtain the spin glass-paramagnetic 

O , phase transition point F/J ~ 1.62 at T = 0. 
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I. INTRODUCTION 

a' 

I The classical random spin systems show various interesting properties that cannot be observed in clean systemSfii^ 

^ The main concern is the existence of the spin glass phase and comprehensive analyses of several models such as the 
Edwards- Anderson'^ and the Sherrington-Kirkpatrick (SK) model^ have revealed the properties of the randomness- 
I t . induced phase transition. 

Ci • However, once we apply a transverse field to the SK model, the model becomes a quantum mechanical one and 
^ , cannot be solved exactly even at the mean field level. Since the work of Bray and Moore^ for the Heisenberg model, 
quantum spin glass models have attracted much interest as they can study the interplay between randomness and 
quantum fluctuations. 

The effect of quantum fluctuations can be easily realized by formulating the models in a path integral form. It is well 
known that the c?-dimensional quantum systems are mapped onto the (d -I- l)-dimensional classical systems.^ In the 
mean-field analysis, models are described by the functional integral of order parameters and the additional coordinate 
("time") dependence demonstrates the fluctuation effects. The simplest possible approximation is to neglect the time 
dependence of order parameters, which misses quantum effects and is not generally justified. 
' Concerning the SK model in a transverse field, the spin glass phase transition has been investigated by many au- 
\^ [ thorsi^'^ i-'^^'-'^^i^^'-'^'^i"'^^'"'^^'^^!"'^'''^^ It was recognized there that the static approximation is not justified at low temperature 
0^ ■ and the spin-glass-paramagnetic-phase transition point at T = is significantly affected by quantum effects. 
CN ' In this paper, we propose a systematic method for quantum random spin systems. It is based on expressing the 
■ partition function using a path integral form. Using the Trotter-Suzuki decompositioni^ we insert the spin coherent 
' state representation as the resolution of unity. The use of coherent states has great advantage for the resulting 
integrals since the spin integral variables are continuous. Furthermore this method can be applied to the arbitrary 
Hamiltonian. The use of the spin coherent states is inevitable for the quantum Heisenberg model and was indeed 
^ , done in Ref.2ft. Although it is not necessary for the SK model, we stress in this paper that our method is convenient 
• '~j ■ for the systematic calculations. We mention, for example, the semiclassical method and the gauge transformation. 
rN ] We apply our formulation to the transverse SK model. We carefully treat the time dependence of the order 
^ I parameters to investigate the spin-glass-paramagnetic-phase transition. We use several field theoretical methods: the 
■ ■ Landau expansion, the renormalization-group-like method at finite temperature, and the derivative expansion at zero 
temperature. The result of the phase diagram is compared to the previous works. 

The paper is organized as follows. In Sec[TT] we introduce a path integral form for the general quantum spin 
Hamiltonian and discuss several methods to calculate the functional integral. In Sec lIIII we apply the method to the 
transverse SK model. After discussing the Landau expansion, we consider improved methods to treat the quantum 
fluctuations. We calculate the spin glass order parameter and determines the phase diagram. Section [IVI is devoted 
to discussions and conclusions. 
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II. FORMULATION FOR QUANTUM SPIN SYSTEMS 
A. Coherent state path integral representation 

The path integral representation for spin systems can be found in, e.g., The closure relation 

is expressed by using the spin coherent states and is inserted to the Trotter-Suzuki decompositioni^ of the matrix 
element of the time evolution operator. The spin operators in the Hamiltonian are replaced by the classical variables, 
and the complex phase factor is added to ensure the dynamics. 

The spin coherent state for spin S is defined by 

|S> =e-^^^3e-*^^=|5), (1) 

where S — (Si, S2, S3) are spin operators satisfying 

[S,,S,]=ie^,kSk, S2= 5(5 + 1), (2) 

and 1 5) denotes the eigenstate of 53 with the eigenvalue 5. 6 and ip are real parameters and parametrize the 
coordinates on a unit sphere. The expectation value of the spin operators is given by 

S = (S I S I S) = 5(sin e cos (p, sin 6 sin ip, cos 9) . (3) 

The closure relation is written as 



/ 



dS|S)(S| = ^^^^ / d{cose) I dip\S){S\ = l. (4) 



We consider the finite temperature partition function for the quantum spin Hamiltonian H — H[S]. We use the 
Trotter decomposition and insert the coherent states defined above. We have the matrix element 



S(t + AT)|e-^^^|S(r)) -exp<^ At 



-(S(r)||-|S(r))-(S(r)|7?|S(r)) 



(5) 



where r is the slice index and represents the imaginary time between and /3 — l/fc^T, and At is the slice width 
which must be taken At — > 0. The first term on the right-hand side is pure imaginary and includes the time derivative. 
The second term gives the Hamiltonian with the operators replaced by the classical variables given by Eq.Q. Thus, 
we arrive at the expression for the finite temperature partition function 

Z = JvScxpi^iJ^ dT$(T) - driJ [S(t)] I , (6) 

where 

$(r) =*(S(r)||-|S(T)) =5(^(T)cos0(r), (7) 

and the periodic boundary condition S(0) — S(/3) is imposed. We note that <i> represents the geometric Berry phase. 
The corresponding term in Eq.([5]) is always imaginary irrespective of the real or imaginary time formulations and 
describes the dynamical motion of S(t) on the sphere. 



B. Calculation methods 



As we mentioned above, the advantage of the coherent state representation is that the spin operators are replaced 
by the continuous classical variables, which is crucial to describe the spin dynamics. Various methods developed in 
the path integral formalism can be applied to the present case as well. Here we discuss two useful methods used in 
the following sections. Without loss of generality in the present context we can confine our discussion to the one-body 
Hamiltonian 

H = -BS, (8) 

where B is the external magnetic field which may depend on time. In the following application of a spin glass model, 
the system can be written in the one-body form by introducing the auxiliary variables. 
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1. Semiclassical method 



First, we mention the semiclassical method which is known as the standard approximation in the path integral for- 
malism. We assume the main contribution comes from the stationary configuration of spin variables. The assumption 
is justified when h = and hence the name semiclassical. The stationary phase equation is nothing but the classical 
equation of motion 

For the Hamiltonian ([5]), it is known that this approximation becomes exact in the coherent state representation. 
The derived path integral representation is ill-defined and we need a regularization to perform the integral explicitly. 
Klauder— used the Wiener regularization and considered the path integral using the stationary phase (saddle point) 
approximation. It was proved, after 20 years, that the approximation becomes exact for the one-body Hamiltoniani^i 
The problem is reduced to solving the classical equation of motion ([9]) under the arbitrary boundary condition. 
The reason why the stationary phase approximation becomes exact is that the coherent states satisfy the minimal 
uncertainty relation. 

In the following application discussed in the next section, this method turns out to be useful when we calculate 
the correlation function. For example, we consider the spin-1/2 system with the Hamiltonian H = — 2r5'3. The 
correlation function in the perpendicular (say, y) direction to the magnetic field T can be written as 

_ , _ JVSny{T)ny{T')exp[iJ^^T)+lrn,{T)] 

where S — Sn — n/2 and — J^^ dr. The result of the integration is given by 

/3r-2r|r-T'| , /3r+2r|T-T'| 
Dir-r') = I ^_,p • (11) 

This form is known and can be obtained by using other methods such as the transfer matrix method. Competitive 
advantage of the present method arises when we generalize the calculation to higher order correlation functions and 
higher spins with S > 1/2. 



2. Gauge transformation 



Second, the gauge transformation is utilized when we calculate the correlation functions. The basis of the state to 
be inserted into the Trotter decomposition can be changed to other arbitrary basis. We consider the rotation in spin 
space by the unitary operator 



Umr)] = exp 



The partition function of the Hamiltonian ([8]) can be written as 



Z = JvSexp i J $(t) + J B(t) • S(t) 



(12) 



(13) 



We consider the time-dependent rotation diagonalizing the Hamiltonian to H = — |B(t)|S'3. Such a choice is always 
possible, but this does not solve the problem because we have the expression 



Z ^ VS exp 



$ (r) + 0(t) • S(t) 



ljBiT)\S,{r) 



(14) 



Due to the time-dependent gauge transformation, the phase acquires an extra term. This cannot be solved generally 
except simple cases such as constant fields and oscillating fields. This type of the gauge transformation was discussed 
in Ref[2^ and the extra term of Eq. (jl4p is called the geometric phase. This expression can be utilized when we consider 
the derivative expansion. When the time-dependence is slow, the geometric phase term is expanded and the only 
remaining thing to do is to calculate the correlation functions as Eq.lll 
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III. APPLICATION TO THE SHERRINGTON-KIRKPATRICK MODEL IN A TRANSVERSE FIELD 

As an application of the spin coherent path integral representation we consider the transverse SK model defined by 
the Hamiltonian with Pauli spins S = cr /2 on lattice sites 

N N 
i^j i=l 

where Jij is the Gaussian random variables with mean Jq/N and variance J'^/N, N is the number of lattice sites, 
and r is the transverse field, i and j run over all points, which means that the interaction is infinite range and the 
mean-field analysis becomes exact. When the transverse field is present, the model cannot be solved exactly and 
we need an approximation. As we explained in the Introduction, the simplest way to proceed is to neglect the time 
dependence of the order parameters introduced as auxiliary integral variables. It looks plausible because the time 
variable is introduced artificially and its dependence cannot directly be observed. However, we discuss in the following 
that the time dependence should be integrated out rather than be neglected, which gives a nontrivial quantum effect. 



A. Mean field theory 



We follow the standard prescription to introduce the auxiliary fields m and q using the Hubbard-Stratonovich 
transformation.^''^ Introducing replicas, we take the average of the nth power of the partition function as 

PSexpK]^ / [^<i>^(r)+^<,(r)] 

^a=l i=l •'^ 

J n N „ t2 " ^ /■ 

+ iT.T. <(-K.W + ^ E E / <Ar)KJr)nUr' 



i<j a, 0—1 i<'j 

The Hubbard-Stratonovich transformation allows us to introduce the order parameters. We have 

72 



X(r') • (16) 



where 



[Z"] = / VmVqc^J^^Yl I K(r)]' - ^ E / [g,fc(T, r')]' + In Tr , 



Tre^ = /pSexp^ / [z<i>''(T)-fr<(r) + Jom„(r)<(r)] + - ^ / g,,(r, T')<(r)n^(r'; 



(17) 



The saddle point equations 

Tr?i2(T)e^ . Trn2(r)n^(r')e^ 

-7^ T , qab[T,T) = 

ire-^ ire 
indicate that m is the magnetization and q the spin glass order parameter 



(18) 



-.(-)= ^^^^c^, ...(r,/)= ^^'"^:i^' (19) 



B. Landau expansion 



To proceed further, we must integrate out the spin variables n"(r) to get the order parameter functional. We do this 
by using the Landau expansion assuming the order parameters are small. The Landau theory for the present model 
was considered in Refl26l by writing down the functional immediately from the symmetry argument. That method is 
phenomenological and the coupling constants in the Landau function cannot be related to the fundamental parameters 
in the original Hamiltonian. Here we derive the Landau function microscopically from the original Hamiltonian. We 
use the derived result to determine the phase boundary where the approximation makes sense. Then the transition 
point can be expressed by T/J and F/J. 

In order to get the result we must calculate the spin correlation function pip . Higher order correlation functions 
can be calculated in the same way as the two point function. For instance the four point function is 



D{ti,T2, rg, T4) = D{t[ -t'^+t'^- O, 



(20) 
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where ti_2.3,4 are arranged in order of magnitude as t[ > T2 > > t^. Higher order functions can be expressed in 
the same way. 

Now we can express each term of the Landau expansion using the correlation functions. For simphcity we consider 
Jo'fTT-aiT) = which means we consider the paramagnetic or spin glass phase. We write for the spin glass order 
parameter 

Qabir, t') = 5abXa{T, t') + (1 - Sab)qab{T, t'), (21) 

to distinguish the role of each term. The first term is the diagonal part in the replica space and represents the spin 
susceptibility. It is unity when F = and the deviation from unity at F ^ represents the quantum effect. The 
second term is the familiar spin glass order parameter. 

Using this representation, we can Landau-expand the averaged free energy as 

/?[/]= -ln(e^^ + e-''-)+limikx:/x^W + ^ E J U^^}, (22) 

where Ik is of fcth order in x and q. For fc = 1 and 2, it is given by 

^1=E/ Xa{Tl,T2)D{n-T2), (23) 

h=^j Xa(Tl,T2)Xa(T3,r4) [D(ti,T2,T3,T4) - D(ti - T2)D(t3 - T4)] 

(J— 1 J T1T2T3T4, 

+ 2^/ qabiTl,T3)qba{T4,T2)D{Tl^T2)D{T3-Ti). (24) 

J T1T2T3T4 

We note that this is the expansion in terms of P'^J^qab{T, r') and (3'^J^Xa{T, r'). Since the order parameters are unity 
at most, we can regard P^J^ as a formal expansion parameter. The order parameters are determined by the saddle 
point equations. In the following we consider the replica symmetric and nonsymmetric cases. 

1. Replica symmetric solution 

We assume that x and q are independent of the replica index. Then the saddle point equations up to first order in 
are given by 

j2 r 

X{ti,T2) = D{ti ~T2) + — / x{t3,T4) [D {ti , T2 , T3 , T4) - D(ti - T2)D{t3 - T4)] , (25) 
q{Tl,T2) =J^f q{T3,T4)D{Tl - T3)D{t4 - T2). (26) 

JT3T4, 

At the linear approximation x(t, r') is equal to the correlation function D(t — t') and the assumption of the replica 
independence is plausible in the perturbative calculation. On the other hand, the time dependence of x cannot be 
neglected. Concerning q, the static approximation seems to be appropriate as we discuss below. 

X can be solved iteratively while q is solved by considering higher order nonlinear terms. The phase boundary can 
be determined by the leading term in Eq. (|26p . Using the static approximation for q, we obtain 

This can be easily solved to find T/ J — 1 a,t T — This perturbative solution is compared to the static approxima- 
tion result F/J = 2 in RefsH and 13 where x is not expanded and treated nonperturbatively. The difference comes 
from the fact that x does not vanish at the phase boundary and contributes to Eg. ((27)) . 
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2. Replica symmetry breaking solution 

It is well known in the classical SK model without transverse field that the spin glass order parameter depends on 
the replica index and the replica symmetry breaking (RSB) solution proposed by Parisi^i is the exact one. Here we 
consider the effect of the transverse field for the RSB solution. The calculation can be done explicitly if we use the 
Landau expansion. The free energy is expanded in qab up to the fourth order and the saddle-point equation is solved 
analytically under the assumption of the RSB.^'^^ This can be done even if the transverse field is incorporated. Near 
the transition point T = Tc, we obtain the expression 

m = Imv Uhy^iQabY ~\C, V QabQbcQca-^C^y^iQabf 

+ \c3j2iQabf{Qacf-lc4Tr{QA, (28) 

abc 

where Qab = P^J^Qab and 9 = [T^ — T^jT^. and the coefficients Ci.2,3,4 depend on x a-nd are expressed perturba- 
tively. We have for Ci^2 

T'lr ^(^i'^2)[i?(ri,r2,T3,r4)-D(ri-T2)D(r3-T4)] + ---, 
C2 = [D{Ti,T2,T3,n)DiT5,T6,T7, Tg) - 3i?(Ti - r2)i?(r3 - n)DiT5 - t^)D{t7 - rg)] + • • • . (29) 

Tc and C3.4 are expressed in a similar way. When F = 0, all the correlation functions are set to unity and we obtain the 
classical result Ci.2,3.4 = 1. Since the effect of the transverse field is only to change the coefficients of the expansion 
(PS)) . the saddle point equation for Qab is easily solved in the same way as the classical case as 

9(2;) = ^? (0 < X < Si = 2I0IC2M), 
02 ^ 

q{x) = \e\ ixi<x< 1), (30) 

where x is the replica continuous variable at n 0. Thus, the transverse field changes the slope of the line. Each 
term of the coefficient Ci^2 can be calculated analytically. If we keep only the leading term in Eg. ([25]) we have 

Ci ^tanh/3rV 1 

2- V-J-LJ 



C2 V /3r 



3 f tanh/3r \ ' 

2 \ pr ) 



1 - 



tanh ISr 



This is larger than unity when PV > 0, which means that the RSB solution approaches the replica symmetric solution 
and we expect that the stability of the replica symmetric solution increases. In the following calculations, we consider 
the replica symmetric solution and our attention is mainly focused on the time dependence of x- 



C. Improved Landau expansion at classical regime 

As we mentioned above the naive perturbative expansion gives the transition point F = J at T = fRefflll) and 
the static approximation for x gives F = 2J at T = 0^^^ Furthermore, several analyses using more sophisticated 
techniques showed that the transition point lies between themJ^ii^ In the following analysis, we reconsider this 
problem using a refined field theoretical method systematically. Our method allows us to obtain the phase structure 
not only at the boundary as was done in previous works but also in the whole space. 

Equation (|25p tells us that the order parameter xi^ij ''2) is approximately equal to the correlation function -D(ti — T2) 
in Eg. pTj) . This function has a slow time dependence at /3F <C 1 and fast at /3F ^ 1. Therefore, when the temperature 
is not so low, the static approximation is expected to be valid. On the other hand, it is not justified at low temperature. 

First we treat the classical regime where the time dependence is not so strong. In this case we can separate the 
time-dependent and -independent parts of XaiT) and ^^^(r, r') as 



Xair) =X + Xq(t), qab{T,T') = q + qab{T, t') , 



(32) 
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where x ^-nd q are zero modes defined as the zero frequency part of the Fourier transformation. Xa(T) and qab{T, t') 
are expected to be small and we consider the Landau expansion in terms of these variables. On the other hand x and 
q are not expanded and treated nonperturbatively as was done in the static calculation. We expand the averaged free 
energy in terms of the nonstatic modes and integrate those modes as 



JvqVxVqVxexp{~Nnf3f[q,x,q,x]} = J VqVxexp{^NnPUs[q,x]} ■ 



(33) 



/off is defined as the classical free energy renormalized by the quantum part of the order parameters. / is expanded in 
X and q up to second order and we carry out the Gaussian integrals. Thus the static and nonstatic parts are treated 
separately and we can take full advantage of both the Landau expansion and the static approximation. 
Using the separation of order parameters, we have 

Tre^ = JvSDzY[Dzaexp ^ / [z^-^Cr) + r<(r) + J(V9^ + V3F^^aX(r)] 

+ vE/ Xair,r'M{r)K(.r') + :^Y. [ 9afc(r, r'K(r)n^(r')l , (34) 



where we introduced the auxiliary variables z and Za (a = 1; 2, . . . , n), and the integration measures are given by 



Dz = 



dz 



Dz 



dz„ 



(35) 



The last two terms are expanded up to second order as 

/" n r t2 " p2 /■ 

Dz\{DzaX{{eP''- + e-P''-)ei,v yT^J;^ drdr' Xa{r,T')DhST ~ t') 
a=l a=l L a=l " •' 

+ / Xq(ti,T2)Xq(t3,T4) (Ti,T-2,r3,T4) - Ai„ (t"i , ^2 ) Ai„ (^3 , T4)] 

a "'-ri,2,3,4 

+ / iab{Tl,T2)qab{Tz,T4)Dh^{Ti,T'i)DhST2,Ti) — --^ , 



where 



K = y^F2 + J2(^z + ^/x~qzaY, 



(36) 



(37) 



and DhST,T') is the correlation function (fTT|) with F replaced by ha- We used the gauge transformation to "diag- 
onalize" the Hamiltonian. In the present case, the magnetic field is independent of time and the extra phase factor 
does not arise here. Then taking the n ^ limit we obtain 



lim^lnTre^ = J Dziln j Dz2{el^^ + e'l^^) 

Y E / Xa(T, r') / Dz,D'z2Dy,{T - t') — 



- lim — 



a •^'^1,2,3,4 J 

+ YE/ gah(ri,r2)gafc(T3,T4) y" Dzx Q D' Z2DhiTi - T3)^^ (^j D'Z2D 



F2\ 

.(r2-r4)^l 



where h is equal to ha with z (z^) replaced by zi (2:2), and 

D'Z2 (• 



„ dZ-\ 2 ^2 /r, 

\/27r 



/ Dz2 (• • Ocosh/?/! 
J Dz2 cosh/3/i 



(38) 



(39) 
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The integrations of the nonzero modes Xa{T) and qab{T) are easily carried out to obtain the renormahzed effective 
free energy with the leading nontrivial contribution 
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Dzi In 



■E 

m— 1 



DziD'z2Dh{m) — 



(40) 



where Dh{rn) is the Fourier transformation of Dhir) and is given by 



tanh [3h 



(/3/i)2 + (7rm)2 (3h 



(41) 



Up to the second order expansion of the nonzero modes, q fluctuations do not contribute to the result. This is because 
the spin glass order parameter is defined as the global variable in terms of the lattice site index as = N f[q\. A 
different conclusion is obtained if we define q as local order parameters qi and write F[q] = f[qi\. Fluctuations in 
each local site give a nontrivial result. Our model was defined as the infinite range model. In that case, q fluctuations 
are not important and the static approximation for q is justified. In this sense the mean field theory of the infinite 
range model is different from that of the finite range model as discussed in Refii. 




0.8 1.2 

r/j 



2.0 



FIG. 1: Phase boundary from the saddle point equations (|42|) . The maximum frequency value of the summation is denoted 
with lines, "m = 0" corresponds to the static approximation result and "Perturbative" denotes the naive Landau expansion 
result (f27)l . "Quantum" is the result in Sec lIII Dl The quantum regime is shown by the shaded area. 



The effective free energy (1401) is a function of q and x- Their values are obtained by solving the saddle point equations 
dfcff/dx = and dfcs/dq = 0. The phase boundary is determined by 9/cff/9xlg=o — and d"^ fcs /d^q\q=o — 0. 
Using formulas derived in the Appendix, we obtain 



(3^J^X 



D'zz 



{PJX? = 



2 / D'zz'^ - 3 ' 



It 



F2 

D'zDnim) — 



D'z {z'~ I D'z'z'^ ) b^im)^— 



(42) 



where h = yT^ + J^x^- 

We show the numerically solved result in FigU] The frequency summation is restricted to a finite value and we 
denote in the figure the maximum value of the summation. The line to = corresponds to the static approximation 
result shown before and "perturbative" means the Landau expansion result ((27)) . We see up to T/J ~ 0.2 from the 
above the maximum value to = 4 is sufficient to find the convergence. The equations cannot be solved at lower 
temperature. Actually we can show analytically that the equations (|1^ do not have a real solution at T = 0. The 
extrapolation of the finite-T results gives F/J ^ 1.54 at T ^ 0, which is close to the results obtained in previous 
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worksjiSii&iil We also show the numerical result of the order parameters in FiglS] The behavior of x shows that the 
quantum effect becomes important at low temperatures and large transverse fields as expected. We also see x ~ 9 at 
low temperature, which implies the reduction of the effective number of the order parameters. 




FIG. 2: Spin glass order parameters q and x from saddle point equations. The frequency summation is taken up to m = 4. 
The bold line denotes the phase boundary. 



D. Improved Landau expansion at quantum regime 



At low temperature the static approximation is not justified and we must use a different method. Introducing the 
auxiliary variables for the expression (j34p . we write the last term in a linearized form 



exp 



t2 » /• 

x{ry)<{T)n:{T') 
\Y.[ Ur)x-\r~T')UT') + jY, j UT)K{r) 



T>z exp 



(43) 



As we analyzed in the Landau expansion x(''") decays exponentially in time at r = and (3. This decay rate is large 
at low temperature and we can use the instantaneous approximation. Using the derivative expansion for z{t), we 
obtain 



Vzexp 



(44) 
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where 



Xi-r), X2 



f3 



(45) 



%2 can be written as X2 = ~c^^x('^)U=o ^tnd we see that this contribution was not taken into account in the previous 
approximation. Assuming the form x(t) ~ e~^^ + e~^^ iP-'^) we can write X2 ~ (/3x)^/2- The auxihary variables 
are integrated out to give 



exp 



exp 



2 



J2 I G(t-t'K(tK(t') 



where 



(46) 



(47) 



Compared to Eg. ([43)1 we find that xi''') is replaced by G{t). Then the correlation function of n1{T)n'^{T') in Eq. ([46l 
is calculated by the gauge transformation as Dii^{t — T')T^/h^. Taking the limit n 0, we obtain approximately 



/3/cff 



Dzi In 



f r2 



(48) 



This effective free energy is valid at low temperature and the phase transition point at T = can be determined from 



pJx 

1 = ■ 



pJx 
1 



I ^^2 



{PJxY 



D'zz 



D'zz^ - 1 



PJ 



D'z z 



D'z'z'^ 



1 



/l2 1+Phx 



(49) 



These equations have the solution (3Jx ^ 0.62 and T/J ^ 1.62 at T — > 0. We see in Fig[T] that the result of the 
quantum regime is smoothly connected to that of the classical regime. 

Finally we determine the boundary between the classical and quantum regime. According to Eq. ()47|) . the time de- 
pendence becomes important when the first Matsubara frequency is comparable to the "mass" term of the propagator. 
Thus we identify the quantum regime ttx ^ 1 which is illustrated in Fig[TJ 



IV. CONCLUSIONS 



We have discussed quantum random spin systems using the field theoretical method based on the spin coherent 
state representation. The partition function is represented as a functional integral and the continuous integration 
variables describe the spin motion on a unit sphere. This formulation can be applied to arbitrary Hamiltonians with 
arbitrary spin S. 

In previous works for the Ising spin {S ~ 1/2) systems, the eigenstates of 5*3 has been used for the closure relation 
to be inserted into the Trotter decomposition. This formulation gives discrete integration variables and the continuum 
limit in the time direction cannot be taken since the time derivative for the discrete variables is ill defined. In this 
sense our formulation is natural and useful even for the Ising systems. 

Our method is also useful for explicit calculations. We can use various field theoretical techniques such as the 
semiclassical method and the gauge transformation. As an application we considered the transverse SK model. The 
classical effective free energy renormalized by quantum fluctuation effects are expressed in terms of order parameters 
and the saddle point equations are solved to obtain the phase diagram. We showed that the time dependence of x is 
important to obtain the result. We found the quantum phase transition point located between the perturbative and 
the static results. Our estimate is F/J = 1.62 at T = and is close to the values obtained by others.— li^i^^ 

What is conceptually important in our calculation is that the role of the order parameter variables is distinguished 
between the static and nonstatic parts. The order parameter is defined as the static part of the variable and the 
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nonstatic part are integrated out to find the effective classical theory, which is reminiscent of the renormalization 
group theory. Therefore, it is a straightforward extension to consider the renormalization group calculation as was 
discussed in Refl26l. Allowing for the spatial fluctuations of the order parameters, we can examine the stability of the 
critical states and calculate the critical exponents. 

We clarified in the present paper the role of the quantum fluctuations using the simple transverse Ising spin glass 
model. It is a straightforward task to apply our results to other models. For example we can consider the transverse 
SK model with arbitrary spin S. In that case, it is not difficult to calculate the correlation function corresponding to 
Eg. ljlOp . When S is large, we find that the time dependence becomes weak and the static approximation becomes a 
good one. In a similar reason, the static approximation is justified when we have an infinite many-body interactioni^ 
These observations show that the present transverse Ising spin model is the simplest one but the quantum effect is 
maximum. We think the next simplest nontrivial application of our method is the quantum Heisenberg model. Most 
of the previous works relied on a semiclassical method such as the static approximation^i^^ and the large- limit^S. 
We hope that our approach will be useful for studying the quantum effects. 

Finally we mention another possible application. In the present paper we considered the imaginary time formulation 
to calculate the partition function. It is also possible to consider the real time formulation which allows us to analyze 
the dynamical correlations. This can be done by using the Keldysh formulation.'^*' We can calculate the dynamical 
correlation function without using the analytic continuation from imaginary to real time. The Keldysh method is also 
useful when we consider the random averaging and field theoretical methods were developed for disordered Fermion 
systems^Si The application to the random quantum spin systems is an interesting problem and is left for future work. 
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APPENDIX: DERIVATION OF THE SADDLE POINT EQUATIONS 



We consider the derivative of the following functions to derive the saddle point equations: 

jDz2fih)gih) 



F = Dzi In 



Dz2f{h) 



G= Dzi 



jDz,f{h) ' 



(A.l) 



1/2 

where h — (V^ + M^) with M = J{^zi + ^/x — 922). After taking the derivative with respect to x or q, we take 



the limit q = 0. First, we consider the derivative with respect to %. We have 



OF 
dx 



Dzi 



= 2 / ^^1 



jDz^fih) 
jDz,^^£-^^f{h) 



2 



Dzi 



jDz2f{h) 

jDz2£.f{h) 
jDz,f{h) ' 



(A.2) 
(A.3) 
(A.4) 



where we referred the definition of the integration measure ([35]) to use ze ^ = — de ^ 1"^ jdz and the partial 
integration in the second line. In the same way, we have 



dG__.V_ 
dx ~ 2 



Dzi 



jDz2£^f{h)g{h) jDz2^f{h)jDz2f{h)g{h) 



jDz2fih) 



jDz2fih) jDz2fih) 



At the limit q ~ 0, h becomes independent of zi and we obtain 



dF 
dx 
dG 
dx 



1 J Dziz"^ -l)f{h) 
" 2^ jDzfih) ■ 



q=0 



1 

2^ 



/ Dziz^ - l)f{h)g{h) J Dz{z^ l)f{h) J Dzf{h)g{h) 



JDzfih) 



JDzfih) JDzfih) 



(A.5) 

(A.6) 
(A.7) 
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where M = J^/\z. We replaced the derivative with respect to M by that with respect to z and used again the partial 
integration. In the same way, we obtain 



dF 
dq 



dG 
dq 



dq^ 



9=0 



9=0 



9=0 



9=0 



1 

1 

'2^ 



jDzzfih) 



jDzfih) 

jDziz^-l)f{h) 
JDzfih) 



0, 



1 jDz{z'-l)f{h) 
JDzfih) 



J Dzz^fih)gih) J Dzz^fih) J Dzfih)gih) 



JDzfih) 



JDzfih) JDzfih) 



(A.8) 

(A.9) 
(A.10) 
(A.11) 
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